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We study a system of few fermions in a one-dimensional harmonic trap, and focus on the case of 
dipolar majority particles in contact with a single impurity. The impurity is used both for quenching 
the system, and for tracking the system evolution after the quench. Employing exact diagonalization, 
we investigate relaxation and thermalization properties. In the absence of dipolar interactions, the 
system is near integrability, and the dynamics remains oscillatory even on long time scales. On the 
other hand, repulsive as well as attractive dipolar interactions lead to quick relaxation to the diagonal 
ensemble average which is significantly different from corresponding thermal averages. A Wigner- 
shaped level spacing distribution indicates level repulsion and thus chaotic dynamical behavior due 
to the presence of dipolar interactions. 


I. INTRODUCTION 


From the microscopic point of view physical systems 
are governed by the laws of quantum mechanics. For sys¬ 
tems of macroscopic size, however, a description on this 
level usually fails due to the huge amount of different de¬ 
grees of freedom. In such case, a tractable approach is to 
take appropriate averages over the degrees of freedom, re¬ 
placing the microscopic laws by the laws of statistical me¬ 
chanics. This switch has some paradoxical consequences, 
such as the emergence of irreversible dynamics from mi¬ 
croscopic laws with time-reversal symmetry. In order to 
reconcile the physical laws on the microscopic and the 
macroscopic level, it seems to be necessary to extend our 
microscopic understanding towards larger system sizes. 
Therefore, we need to improve our capabilities in obtain¬ 
ing exact or quasi-exact knowledge about a system. 

Maybe the most promising direction to achieve this 
goal are quantum simulations, that is, experiments with 
well-controlled, tunable quantum systems lj. A re¬ 
markable step towards bridging the gap between micro¬ 
scopic and macroscopic physics were the recent experi¬ 
ments in Heidelberg with cold fermionic atoms in a one¬ 
dimensional trap . Outstanding through the precise 
control over particle numbers, the experiment of Ref. [4j 
demonstrated the formation of a Fermi “sea” of up to five 
majority atoms interacting with a single impurity atom. 
These experimental advances have also boosted the the¬ 
oretical interest in one-dimensionally trapped Fermi sys¬ 
tems @-[l3|- 

In this article, we extend the scenario from the Hei¬ 
delberg experiments into two directions: First, we as¬ 
sume dipolar interactions between the fermions. Sec¬ 
ond, we shift the focus onto dynamical processes. The 
“standard” interactions in quantum gases are contact in¬ 
teractions. With this, a spin-polarized Fermi sea as in 
Ref. [H remains non-interacting, since contact interac¬ 
tions are forbidden by the Pauli principle. Clearly, the 
most prominent fermions in nature are electrons with 
long-range Coulomb interactions. To design an atomic 
quantum simulation of a Fermi systems with long-range 
repulsion, one might stick to dipolar interactions. The 
long-range character of these interactions is certainly one 


of various reasons for the huge research interest in dipo¬ 
lar quantum gases 1.8-26j. With advances in cooling 
molecular systems EzTI l 28 j, and recent developments in¬ 


volving Rydberg gases 


new experimental plat¬ 


forms involving strong dipolar interactions are emerging. 
One-dimensional systems with dipolar interactions can 
provide a particularly interesting scenario, as the confine¬ 
ment modifies the effective interaction potential with the 
possibility of confinement-induced resonances 


In a system which is strongly confined to one dimen¬ 
sion, dipolar interactions become almost local, and there¬ 
fore barely modify the ground state of a small Fermi 
gas. Nevertheless, the dynamical behavior of the sys¬ 
tem might be affected strongly. Accordingly, our paper 
is concerned with the dynamics of the system. Cer¬ 
tainly, from the point of view of numerics, dynamics 
is the most difficult aspect of a quantum system, and 
therefore best-suited to be studied via a quantum sim¬ 
ulation. Quantum systems out-of-equilibrium, and the 
relation between unitary time evolution and thermody¬ 
namics, has attracted broad research interest, cf. Refs. 
fTil - l42l or the review articles [43!, [ 3 . A closely related 
subject are Lieb-Robinson bounds and the velocity 
of correlation spreading, in particular in the presence of 
long-range interactions I16H~)2I |. 

In this paper, we will study a one-dimensional Fermi 
system in which majority atoms interact, via contact, 
with a single impurity, as in the Heidelberg experiments. 
Among themselves, the majority atoms can interact via 
dipolar interactions. We will first provide a brief descrip¬ 
tion of this system, and introduce the effective dipole- 
dipole interaction in Section |H] In Section IIIII we will 
sketch the concept of time evolution in quantum me¬ 
chanics, and define some useful quantities to describe the 
dynamics of a many-body system. After briefly studying 
static properties of the ground state in Section[[Vl we will 
focus on relaxation processes from an out-of-equilibrium 
state in Section El Here, the impurity will serve, in a 
two-fold way, as a control knob: On the one hand, we 
can excite the system by exciting the impurity atom to 
higher oscillator levels. On the other hand, the impurity 
can serve as a “thermometer”. We will investigate under 
which conditions the energy of the impurity atom relaxes 
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to a constant value, and whether this value corresponds 
to the thermalized value. It will turn out that dipolar 
interactions between the majority atoms significantly en¬ 
hance relaxation. Generally, the diagonal ensemble de¬ 
scribes well the relaxed state, but it differs from thermal¬ 
ized values. Finally, in Section fVTl we will relate our re¬ 
sults to more general properties of the system and analyze 
the level spacing distribution. Although not an integrable 
model, the system without dipolar interactions exhibits 
a strong tendency of level clustering. On the other hand, 
a Wigner-sliaped distribution for the system with dipolar 
interactions indicates level repulsion and chaotic dynam¬ 
ics. In this sense, the level spacing distribution turns out 
to be a good indicator for the system’s relaxation proper¬ 
ties. In the appendix, we calculate the interaction matrix 
elements of a dipolar gas in one dimension. 


II. SYSTEM 

We consider a system of fermionic atoms in an effec¬ 
tively one-dimensional harmonic trap Vt ra p = ^-^x 2 , 
where m is the mass of the atoms, and u> the axial trap¬ 
ping frequency, chosen along the ^-direction. The one- 
dimensional trapping is achieved via a sufficiently strong 
anisotropy of the three-dimensional trap, freezing out the 
dynamics in the transverse directions. The ratio between 
axial trapping frequency to and the transverse trapping 
frequency uj± is given by the parameter A = y/u/ui±_ = 
l±/l , with l = ■) and l± = y/h/(muj±) the cor¬ 
responding harmonic oscillator lengths. For most of our 
numerical studies, we choose A = 1/3, the value reported 
in Refs. 01- We will work in harmonic oscillator units 
of the axial trapping, that is,w = l,Z=l,m=l, and 
h= 1. 

The fermions can be prepared in different internal 
states, or, equivalently, the system may consist of differ¬ 
ent species with equal mass. We restrict ourselves to two- 
component systems, where A/- and denotes the num¬ 
ber of fermions in each component. We focus on the case 
where several majority atoms (forming the t-component) 
interact with a single impurity atom (N± = 1) via contact 

interactions, Kontact = 9 (Eieto'et ^ Xi ~ x i) witl1 effec “ 
tive interaction strength g. Such scenario has been stud¬ 
ied experimentally in Ref. 0. Since identical fermions 
do not interact with each other via contact interaction, 
the Fermi sea of the majority systems would be non¬ 
interacting. To overcome this, we may equip the major¬ 
ity atoms with an electric or magnetic dipole moment 
d, which leads to dipole-dipole interactions between the 
majority atoms. For simplicity, we assume that the impu¬ 
rity atom carries no dipole moment, such that it interacts 
with the majority atoms only via contact interaction. 

The dipoles shall be aligned, forming an angle 9 with 
the cc-axis. Then, the scattering between dipoles, sepa¬ 
rated by a vector r, is described by the usual interaction 
potential Vdd(r) = 73(1 — 3 cos 2 9 r d), with 9 r d the angle 
between the dipole and r. Due to the strong transverse 


confinement, one can assume that transverse degrees of 
freedom will not be excited, and integrate out the y an d 
z variables. This leads to the effective potential [321-1 34 1 

T/eff/ n d 2 [l + 3 cos(20)] 

Vdd( x ) = -- U(x/X) = uU(x/X), (1) 

with the dimensionless potential 

U(x) = — 2\x\ + V2 tt( 1 + x 2 )e x ^ 2 erfc(a;/v / 2), (2) 

where erfc is the complementary error function. An ad¬ 
ditional term oc S(x) occurring in the effective potential 
can be neglected for identical fermions. We note that the 
strength of the dipolar interactions, u, can be tuned via 
the polarization angle 9 , and can be made attractive or 
repulsive. 

The total Hamiltonian will be studied by exact diag- 
onalization in the Fock basis of the harmonic oscillator 
levels. We truncate this basis, taking into account up to 
22 single-particle levels, for systems of up to five majority 
atoms. We consider interaction strengths up to g = 10 
and |u| = 10. In second quantization, the Hamiltonian 
reads (assuming no dipole moment for the impurity j. 
atoms): 

H=J2 + f E 

+ 5 E Kir( 3 ) 

"" (jjfcz)et 


with annihilation/creation operator aj and aj, where the 
index j refers to the jth orbital, while the index for 
the spin has been suppressed. A constant energy term 
Nhuj/2 , corresponding to the zero-point energy of the 
harmonic oscillator with N particles, is neglected. With 
ifij (x) denoting the normalized orbitals, the interaction 
matrix elements are defined as follows: 


V r contact 
ijkl 


and 


t rdipolar 
V ijkl 


dx dx'tpi(x)(pj(x')ipk(x)ipi(x') 



dx dx' XJ{x/X 


x'/X) 


X ipi(x)ipj(x')ipk(x)ipi(x'). 


( 4 ) 


( 5 ) 


Both integrals can be solved by integrating out the 
center-of-mass coordinate x + x'. While the remaining 
integral in the relative coordinate x — x’ is trivial for the 
contact interaction, some care must be taken for the dipo¬ 
lar potential. We describe how to evaluate this integral 
in the appendix. 


III. TIME EVOLUTION AND 
THERMODYNAMICS 

In this section, we will define different quantities which 
are useful for analyzing the dynamical properties of a 
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quantum system. We will be interested in scenarios 
where the system is initially prepared in some pure state 
|^ini), and then evolves in the Schrodinger picture under 
the time-independent Hamiltonian H. The time-evolved 
state is then given by 

Mt)) = z~ lHt l^ini) = E |a). (6) 

Ot 

Here, we have decomposed the initial state into eigen¬ 
states |a) of H , H |a) = E a |a). The coefficients c a are 
given by c a = (a| H'mi). 

Using exact diagonalization, we can determine the full 
energy spectrum of a small system, and then straight¬ 
forwardly calculate the evolution of an arbitrary state. 
The state vector itself, however, is not accessible ex¬ 
perimentally. In practice, the system dynamics can 
be tracked by looking at the evolution of any observ¬ 
able. From the computational point of view, this re¬ 
quires to evaluate the quantum average of operators O: 
(0(t)) = <tf(t)|0|tf(t)>. If all eigenstates a are non¬ 
degenerate, ( 0(t )) will, for sufficiently long times, con¬ 
verge to a value described by the diagonal ensemble, Ode: 

(0{t)) =J2 c 0 c «e- i[E '*- E e )t (P\ O |a) 

a,/3 

—> ^ \c a \ 2 (a| O |a) = Ode- (7) 

a 

In the limit on the right-hand side, any off-diagonal con¬ 
tributions (a | O |/9) have been averaged away by the os¬ 
cillatory factor. 

In contrast to the diagonal ensemble, which may keep 
some memory of the initial state, thermal ensembles de¬ 
fine the state of a system only by a few thermodynamic 
parameters which are fixed by the initial conditions. For 
example, if the canonical ensemble is applied, the state 
of the system is defined by the temperature T. On the 
other hand, to associate a temperature with the state of 
the system we shall evaluate the ensemble average E t h(T) 
of the energy: 

£th(T) = ±Tr , (8) 

where the Boltzmann constant has been set to 1, and Z 
denotes the partition function Z = Tie~ H / T . Compar¬ 
ison of the thermal energy E^T) with the system en¬ 
ergy, E = (’Fi n i| H I'Fini), yields an effective temperature. 
With this, one is able to calculate thermal expectation 
values of any observables: 

Oc E (T) = iTr(Oe- ff / T ). (9) 

The observable which, in the present study, will be 
used to track the system is the harmonic oscillator level 
of the impurity. Experimentally, this quantity can be 
determined by reducing the height of the trap. At a 
certain height, the impurity escapes the trap, and from 


this value, the oscillator level of the impurity can be de¬ 
termined. Formally, the corresponding operator reads 
m imp = Y /m ’m\n rn ) l (n m \ l . It acts only on the impu¬ 
rity (that is the .(.-component), and counts the number 
of particles n m in each orbital m , weighted by the level 
number. We can associate (mi mp ) with the energy of the 
impurity. 


IV. STATIC PROPERTIES 

Before studying the relaxation dynamics of the system 
in the next section, let us first discuss some static proper¬ 
ties. As has been demonstrated in Ref. [H, it is possible 
to measure the interaction energy of the impurity with 
high accuracy. The measurement is based on determin¬ 
ing the resonance frequency for an internal transition of 
the impurity atom, which, in the presence of majority 
atoms, is shifted by the amount of the interaction en¬ 
ergy. Increasing the number of majority atoms from one 
to five, it was shown that the interaction energy quickly 
approaches the theoretical predictions made for an im¬ 
purity in a homogeneous Fermi sea p53| . Here, we will 
equip the majority atoms with a dipole moment, that is, 
the Fermi “sea” becomes interacting. As we show below, 
the interaction energy of the impurity provides informa¬ 
tion also about the strength of the dipolar interactions, 
although they are restricted to the majority atoms. 

The system is described by the Hamiltonian H of Eq. 
©• Two control parameters, u and g , allow to ad¬ 
just independently the interaction between the majority 
atoms, and the interactions with the impurity. Using 
numerical diagonalization we have determined the total 
energy E(u,g) for different dipolar and contact inter¬ 
action strengths. We obtain the interaction energy of 
the impurity with the majority atoms as the difference 
E(u,g) — E(u , 0). As shown in Fig. [T| repulsive interac¬ 
tions between the majority atoms (u > 0) lead to a re¬ 
duced interaction energy of the impurity, while attractive 
dipolar interactions (u < 0) will increase the impurity’s 
interaction energy. 

An explanation for this behavior are density changes 
due to the dipolar interactions: Repulsive interactions 
reduce the local density of majority atoms near the trap 
center. Clearly, this also leads to decreased spatial over¬ 
lap between majority atoms and impurity, since the im¬ 
purity predominantly occupies the lowest oscillator level, 
and thus has an amplitude maximum at the center. As 
a consequence, the contact interaction energy between 
both species is reduced. Oppositely, attractive dipolar in¬ 
teractions lead to an increased density of majority atoms 
in the center of the trap, and therefore enhance the con¬ 
tact interactions. 
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FIG. 1: (Color online) Interaction energy of an impurity 
atom with several majority atoms as a function of the sys¬ 
tem size. The black line shows the behavior for a system 
of non-interacting majority atoms, u = 0. Repulsive dipo¬ 
lar interactions between the majority atoms, u > 0, increase 
the interaction energy, whereas attractive dipolar interactions 
decrease the interaction energy. 


V. DYNAMIC PROPERTIES - CASE STUDY 

We will now investigate how the dipolar interactions 
between the majority atoms modify the properties of the 
system dynamics. Again, the impurity atom provides 
a tool for controlling and measuring the system. More 
concretely, we propose to study the time evolution of a 
system out of equilibrium by tracking the energy of the 
impurity. Numerically, determining the time evolution is 
an extremely hard task, as knowledge of all eigenstates 
is required, while exact diagonalization algorithms work 
most efficiently if only few eigenvectors are obtained. We 
therefore restrict the system size to 3+1 atoms. Taking 
into account parity symmetry of the Hamiltonian, and 
restricting the single-particle basis to 22 (20) states, we 
have Hilbert space blocks of 16940 (11400) states. We 
have checked our results for convergence with respect to 
increases in the single-particle basis. Clearly, truncation 
errors are typically more pronounced for excited states 
than for the ground state, and time evolution exponenti¬ 
ates every error. Accordingly, in some cases, convergence 
can be achieved only on short time scales. 

We choose an easy-to-prepare initial state at relatively 
low energy. A convenient choice, taken for Fig. (2j is a 
Fock state with the three majority atoms populating the 
three lowest levels. In the absence of an impurity, this is 
the ground state of a system without dipolar interactions. 
The impurity is used to bring the system out of equilib¬ 
rium. To have pronounced effects, the impurity shall 
occupy a level significantly above the equilibrium value. 
We have chosen m i mp = 8 in Fig. [2] We then studied the 
time evolution of this initial state under the Hamiltonian 
H of Eq. © for a contact interaction strength g = 1, and 
for different values u for the dipolar interactions between 
the majority atoms. Being a Fock state, the initial state 
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FIG. 2: (Color online) Time evolution tracked by the average 
level number (m)i mp of the impurity. The initial state is a 
Fock state with the three majority atoms in the three low¬ 
est levels, and the impurity occupying the 8th excited level, 
riimp = 8. In all plots, the Hamiltonian H contains a repulsive 
contact interaction of the majority atoms with the impurity 
(g = 1). The strength u of the dipolar interactions between 
the majority atoms takes the value denoted in each plot. The 
horizontal dashed lines show the corresponding diagonal en¬ 
semble average. The red dotted line indicates whether the 
numerical treatment is converged: The black line shows the 
results taking into account 20 (a-c) or 22 (d,e) harmonic os¬ 
cillator levels. The red dotted shows the results for a smaller 
Hilbert space, taking into account 18 (a-c) or 20 (d,e) levels. 


is an (excited) eigenstate of the non-interacting Hamil¬ 
tonian, g = u = 0. This allows to view the scenario also 
as an interaction quench, during which both contact and 
dipolar interactions are suddenly switched on. 

Even without quantitative analysis, it is obvious from 
Fig. [ 2 ] that relaxation does not take place in the absence 
of dipolar interactions [Fig. ©a)]. In this case, the sys¬ 
tem evolution is characterized by an oscillatory behavior 
for all times. In particular, even after a long relaxation 
time, the amplitude 5(m i mp ) of these oscillations still at¬ 
tains large values, 6(mi mp ) > 1. On average for long 
times (350 < t < 1000), the impurity level ( m\ mp ){t ) 
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FIG. 3: (Color online) Time evolution as in Fig. 0 but with 
strong contact interactions, g = 10, and no dipolar interac¬ 
tions, u = 10. The black line (red dotted line) shows the 
results taking into account 22 (20) harmonic oscillator levels. 
The horizontal dashed line indicates the diagonal ensemble 
average. 

oscillates around a mean g = 3.7, close to the diago¬ 
nal ensemble average of 3.6, with a standard deviation 
cr = 0.9. 

One might argue that the observed regularity of the dy¬ 
namics is due to the weakness of interactions, and there¬ 
fore the proximity of the initial state to an eigenstate of 
the Hamiltonian H. As a quantitative measure for the 
amount of Hamiltonian eigenstates participating in the 
initial state, we define a kind of entropy 

S = |ca| 2 ln|c a | 2 , (10) 

a 

with the c a the coefficients of a decomposition of the ini¬ 
tial state in the eigenbasis of H , cf. Eq. ( 0 . If S = 0, 
the initial state is an eigenstate of H. If S = In D, D 
eigenstates contribute with equal weights. Indeed, for 
the Hamiltonian of Fig. [2] (a), the entropy of the initial 
state takes only a moderate value, S = 2.7. Also, the 
average interaction energy (V) = 0.2 is small compared 
to the total energy, (H) = 11.2. However, even for an 
increased contact interaction g = 10 (and u = 0), the dy¬ 
namics remains regular, see Fig. [3l although now entropy 
and interaction energy are significantly larger, S' = 4.1 
and (V) = 2.4. This suggests that it is not the weak¬ 
ness of the quench which led to the regular dynamics in 
Fig. [2] (a). We note that for g = 10, the time evolution 
is quantitatively not converged for large times t > 100. 
Qualitatively, however, the evolution remains regular on 
all time scales, independently from the number of single¬ 
particle states taken into account. 

Next, we consider the case of weak dipolar interac¬ 
tions between the majority atoms, in addition to weak 
contact interactions with the impurity, see Fig. [2](b,c). 
The oscillation amplitudes are reduced compared to the 
case without dipolar interactions, with a long-time mean 
g = 3.5 ± 0.7 for attractive interactions (u = —1), and 
g = 3.6 ± 0.5 for repulsive interactions (u = 1), in agree¬ 
ment with the diagonal ensemble average (3.5 in both 
cases). However, in particular for the attractive system, 
periods of strong oscillations occur repeatedly even for 



g ± a 

S 

(Vdd) 

g = i,w = o 

3.7 ±0.9 

2.7 

0 

g — l,u — —1 

3.5 ±0.7 

2.8 

-0.4 

g =l,u=l 

3.6 ±0.5 

2.8 

0.4 

g = l,u = —5 

6.0 ±0.3 

2.4 

-2.1 

g = 1, u = 10 

3.4 ±0.4 

4.1 

4.1 


TABLE I: Quantitative analysis of the data presented in Fig. 

m 

long relaxation times, keeping the system away from a 
steady-state. In both the attractive and the repulsive 
case, the entropy of the initial state is S = 2.8. Addi¬ 
tionally to the contact interaction energy, (14) = 0.2, the 
total interaction energy now contains also a dipolar con¬ 
tribution: (V) = (Vdd) + (Vc). The dipolar contribution 
takes the values (Vdd) = ±0.4 for u = =pl- 

In Fig. G](d,e), we turn our attention to strong dipo¬ 
lar interactions, u = — 5 and u = 10. Since attractive 
interactions have a stronger tendency to populate higher 
oscillator levels, we are not able to achieve numerical re¬ 
sults which remain convergent for long times. As seen in 
Fig- [2](d), the time evolution for u = —5 is quantitatively 
not converged for t > 50. This means that only the initial 
decay of (n; mp )(t) is captured accurately. However, for 
larger times the evolution consists only of small-valued 
fluctuations around a mean value corresponding to the 
diagonal ensemble average (g = 6.0 ± 0.3 for a basis with 
22 states, and g = 6.2 ± 0.3 for a basis with 20 states). 
This suggests that the system quickly reaches a steady 
state, despite the remarkably small entropy of the initial 
state, S = 2.4. The interaction energy is (Vdd) = —2.1. 

For strongly repulsive dipolar interactions, shown in 
Fig. [1(e), we are able to obtain fairly well converged re¬ 
sults even for long relaxation times. The system quickly 
evolves towards the average of the diagonal ensemble 
(3.4), around which it fluctuates in an erratic manner 
with a standard deviation a = 0.4. This means that 
the deviations from the average are only slightly smaller 
than for weak repulsive interaction, Fig. [2] (c). Instead, 
it is the absence of oscillatory behavior which make the 
strongly interacting system appear significantly more re¬ 
laxed than the weakly interacting one. Both, average 
dipolar interaction energy and entropy of the initial state 
are large: (Vdd) = 4.1 and S = 4.4. 

Finally, we may ask whether this relaxed state is also 
thermalized. In the canonical ensemble, energy E is as¬ 
sociated with temperature T according to Eq. ©. This 
leads to a gauge curve E(T) plotted in Fig. 0] (a) for 
<7=1, and u = 0 and 10. From this curve, we infer an 
effective temperature T = 3.2 for the system with strong 
dipolar interactions. On the other hand, using the canon¬ 
ical ensemble Eq. ©, we can also evaluate the thermal 
average of rii mp for any given temperature. The question 
is whether rii mp relaxes to the value which corresponds 
to the effective temperature defined via the energy. The 
temperature dependence of rii mp , which has basically no 












6 



0 2 4 6 8 10 

T 


FIG. 4: (Color online) Canonical ensemble averages of (a) en¬ 
ergy Eth and (b) impurity oscillator level (m; mp ), as a function 
of temperature T. 

dependence on u, is plotted in Fig. |U(b). For T = 3.2, we 
find the thermal average (ni mp ) t h ~ 2.8, which is signifi¬ 
cantly below the temporal average in Fig. [2] (e), around 
3.5. This comparison shows that the system, although it 
relaxes, does not thermalize. 

VI. DYNAMIC PROPERTIES - DISCUSSION 

Let us summarize the main findings of the case studies 
presented in the previous section: 

• Relaxation towards the diagonal ensemble average 
occurs in the presence of sufficiently strong dipolar 
interactions between the majority atoms. 

• The system dynamics remains oscillatory for a sys¬ 
tem with purely contact interactions. This seems 
not to be related to a smaller interaction energy, 
nor to a smaller participation entropy. 

• The system, even when it relaxes, does not ther¬ 
malize. 

Let us now try to understand the differences in the re¬ 
laxation behavior from a more general point of view. A 
concept which will prove useful is the distinction between 
chaotic and regular dynamics: While regular dynamics 
preserves the coherences of the initial state, chaotic evo¬ 
lution should quickly lead to relaxation. Then, the differ¬ 
ent behavior should be reflected by the energy spectra. In 
quantum mechanics, chaotic behavior is characterized by 
level repulsion. For a system with time-reversal symme¬ 
try, the Flamiltonian matrix is real-symmetric, and ran¬ 
dom matrix theory, cf. Ref. (54|, tells us that the level 
spacing should follow a Wigner distribution function, 
-Pwigner(s) = \se~ s n / 4 . Clearly, the maximum of this 
distribution is located at finite spacings s = yj2/-K. On 
the other hand, such level repulsion is not expected for 
a system with regular dynamics. In this case, conserved 
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FIG. 5: (Color online) Unfolded level spacing in energy spec¬ 
tra of H with <? = 1 and for different u. For normalizing the 
distribution, we did not take into account levels at s = 0. 


quantities will allow energy levels to cross, that is, corre¬ 
lations between the levels are absent. This may lead to a 
Poissonian level spacing distribution, Pp 0 i S son(s) = e -s , 
with a maximum for s = 0. 

Before analyzing the level spacing distribution, it is 
necessary to unfold the energy spectra. Naturally, level 
spacings will be larger in spectral regions with a low den¬ 
sity of states than in dense regions, but this obviously has 
nothing to do with correlations between energy levels. 
By unfolding the energy spectrum, we shall guarantee 
an overall homogeneous density of states, such that level 
correlations become the only source for variations in the 






















































7 


level spacing. To that aim, we follow the unfolding pro¬ 
cedure described in Ref. j54[ . The first step is to smooth 
the discrete density of states via a convolution with a 
Gaussian: 


i N i 

p sm oME)- ( 11 ) 

z— 1 V 

We have chosen the width A of the Gaussian to be given 
by five times the mean level spacing. The smoothened 
density of states also leads to a smooth staircase func¬ 
tion ^smooth (E), that is, a smooth cumulated density of 
states: 


/ E 

dE' p S mooth{E') (12) 

-oo 

The unfolded energy levels e,; are then associated with 
ei = NT, smooth (Ei) , and the unfolded level spacings are 

Si — 6-i Ci—i. 

We have evaluated the unfolded level spacing distri¬ 
bution for the energy spectrum of H with g = 1 and 
u = 0,10,-10, for 3+1 atoms, taking into account 20 
harmonic oscillator levels. In the presence of dipolar in¬ 
teractions, as seen in Fig. [5] (b,c), the levels spacing is 
well described by the Wigner distribution, except for a 
somewhat increased number of levels near s = 0 in the 
case of repulsive interactions, see Fig. [5] (b). Never¬ 
theless, the distributions clearly indicate level repulsion, 
and thus chaotic dynamics. This is in agreement with 
the relaxation properties discussed in the previous sec¬ 
tion. More complicated to interpret is the distribution in 
the absence of dipolar interactions, Fig. [5] (a). Neither 
the Wigner nor the Poisson distribution capture the level 
spacing distribution. On the other hand, an eye-catching 
property of the distribution is the exorbitantly large num¬ 
ber of levels near s ~ 0, with around 4400 of 11399 spac¬ 
ings peaked around zero. This number is not compatible 
even with a Poisson distribution, for which around 10% 
of the levels should be located within s < 0.1 (57]. The 
distribution of the remaining levels is neither described 
by a Poisson nor by a Wigner function, even if we exclude 
the s = 0 levels from the normalization (which leads to 
a normalized distribution of the remaining levels). The 
Wigner distribution catches well only the behavior at in¬ 
termediate spacing, 1 < s < 2, whereas the tail of the 
distribution has Poissonian shape. It might be that the 
huge number of degenerate levels spoils our unfolding 
procedure, and/or that a mixed phase space, consisting 
of both chaotic and regular regimes, leads to such un¬ 
conventional distribution function. As a conclusion, the 
dominant feature of the distribution are the abundantly 
many degeneracies, and therefore we may take the ob¬ 
served level spacing as an indicator for a regular system 
dynamics. 

An observation which might, at least partly, ex¬ 
plain the large number of nearby levels is the “quasi”- 
integrability of the model without dipolar interactions: 
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FIG. 6: We plot the average occupation rriimp vs. the energy 
for all eigenstates of H , with g = 1 and for different u. 


First, there is an exact solution for 1+1 particles in a har¬ 
monic trap [ssl ]. and a quasi-exact ansatz for N +1 parti¬ 
cle [l6| . Second, in the homogeneous case with N+l par¬ 
ticles, the system is integrable via Bethe ansatz [53L i56l|. 
Finally, although no exact ground state solution is known 
for the trapped system with N+l fermions, exact wave 
functions are easily obtained for several excited states 
due to the so-called “fermionization” of the problem. In 
fact, the degeneracies in the energy spectra can be traced 
back to these fermionized solutions. The idea behind the 
fermionization is to fill N +1 different harmonic oscillator 
levels in a fermionic way, that is, by assuming the wave 
function to be given by the corresponding Slater determi¬ 
nant. Such spatial wave function is fully antisymmetric, 
irrespective of the particle spin. Obviously, it completely 
suppresses the contact interactions, and in the absence 
of dipolar interactions, it provides an eigenstate with the 
energy given by the single-particle energy of the occupied 
harmonic oscillator levels. For larger energies, there is a 
rapidly growing number of ways how N+l particles can 
be distributed, leading to a huge number of degenerate 
levels. In this context, it is interesting to note that for 
strongly repulsive contact interactions, the energies of all 
eigenstates approach the values given by these fermion¬ 
ized solutions (til.[7!. In that case, energy levels are either 
degenerate, or separated by the integer spacing Hco, which 
apparently would lead to a very untypical bimodal level 
spacing distribution. 

Let us finally study why the system, even for cases in 
which it equilibrates, lacks thermalization. A hypothesis 
explaining thermalization of an isolated quantum system 
is the eigenstate thermalization hypothesis. It assumes 
that, for any eigenstate, the quantum average of an op¬ 
erator and the corresponding eigenenergy are correlated, 
that is, all eigenstates in the vicinity of some energy E 
have the same expectation values with respect to rele¬ 
vant observables O. Moreover, these values are assumed 
to be “thermalized”, that is they shall coincide with the 
thermodynamic average. We test the eigenstate ther- 







malization hypothesis by plotting the average occupation 
(wiimp) versus the energy for all eigenstates in Fig. [6] In 
both cases, with or without dipolar interactions, the val¬ 
ues of (rriimp) are spread over a broad range for almost 
any energy in the spectrum. Accordingly, the eigenstate 
thermalization hypothesis does not hold, and fixing the 
energy in a microcanonical ensemble will not fix the ex¬ 
pectation value of to i m p. This explains why thermaliza¬ 
tion does not take place. 

VII. SUMMARY AND OUTLOOK 

In summary, we have suggested to probe the dynamics 
of a Fermi system trapped in ID by exciting an impu¬ 
rity and tracking the oscillator level of the impurity. We 
have studied a system with dipolar interactions, and have 
contrasted its dynamics to the dynamics of a system with 
purely contact interactions. While the latter shows quan¬ 
tum collapse and revival effects even after long relaxation 
times, the system with sufficiently strong dipolar inter¬ 
actions relaxes quickly to its diagonal ensemble average. 
We relate this finding to the level spacing distribution 
which indicates chaotic behavior in the presence of dipo¬ 
lar interactions. In the case without dipolar interactions, 
a huge number of degenerate levels characterizes the level 
spacing distribution. 

Even despite the small system size of only four 
fermions, we encounter cases where the time evolution 
cannot be computed faithfully on long time scales. This 
demonstrates the need for better computational tech¬ 
niques, and/or alternative approaches such as quantum 
simulations. Our work suggests to explore dynamical as¬ 
pects in the experimental setting of Ref. 4]. In view 
of the recent experimental progress with dipolar atoms 
HE HE HU an d molecules HU , scenarios as the one stud¬ 
ied in this paper are becoming experimentally feasible. 
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Appendix A: Dipole matrix elements 

We need to evaluate matrix elements Vyy° lar = 
f^° 0 dx dx' U(x/X — x'/X)ipi(x)ipj(x')ipk{x)ipi(x') for a 
two-body potential 

U(x) = — 2\x\ + v / 27t(1 + x 2 )e x / 2 erfc(a;/v / 2). (Al) 
The orbitals ip, (x) are the usual harmonic oscillator levels 

= \l~rj^jji e ~ x2/2Hn ^ = Nje- x2/2 H n (x), 

(A2) 

where H n (x) are the Hermite polynomials, and Nj de¬ 
notes the normalization constant of the wave function. 
To evaluate the integral, we replace the coordinates x 
and x' by relative coordinates r = x — x' and center- 
of-mass coordinates R = x + x'. The integral over the 
center-of-mass part, I^\ reduces to 

4° = J R i ~ 1 e~ R2 ^ 2 dR = 2*/ 2-1 [1 + (-l) i+1 ] F ^ , 

(A3) 

with i some positive integer. Also the integral over the 

(i) 

relative coordinate, // ; , turns out to have a compact 
solution 




r i ~ 1 e -r / 2 


-2 


V2tt(1 + 7 T?)e^ / 2 erfc 
X z 
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x\ 

V2X 

)\ 


dr = i 2~i-i [-1 + (-1) 1 ] 

A 


2 i+1 r 


i +1 


— 2 v / 7rA i+1 r(i) 2 -Fi ( 2 


i i -hi i + 2 


l-A 2 - v^A l+1 r(i + 2) 2 Fi -, 


i i -hi i 2 


l-A- 

(A4) 


where 2 ^ 1 ( 0 , 6, c; z) is the regularized hypergeometrical func tion. 


Any matrix element V^° lal can be decomposed into with prefactors stemming from the Hermite polynomials 
a sum over products of these two integrals J 1 ^ and 1 ^, 
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and the normalization factors Nj of the orbitals. How¬ 
ever, care must be taken when numerically evaluating the 
matrix element in high orbitals. Since the numeric val¬ 
ues of the integrals become large, while the normalization 
factors become small, it is crucial to keep symbolic ex- 

_I 

integral[nl_, n2_, n3_, n4_,Lambda_] := Module[{Listl, List2, IR, Ir, r, R}, 

Listl = CoefficientList[HermiteH[nl,(r+R)/2]HermiteH[n2,(R-r)/2] 

HermiteH[n3,(r+R)/2]HermiteH[n4,(R-r)/2],R]; 

IR = Table[I_R[i], {i, 1, Length[Listl]}]; 

List2 = CoefficientList [Dot [Listl, IR], r]; 

Ir = Table[I_r[i,Lambda], {i, 1, Length[Listl]}]; 

N[Dot[List2, Ir] norm[nl] norm[n2] norm[n3] norm[n4]]]] 

where I_R[i] is given by /r(z), I_r [i, lambda] is given by I r (i) for some choice of A, and norm [nl] is the normalization 
factor N n \. 


pressions as long as possible to avoid numeric errors. A 
piece of Mathematica code which performs the decompo¬ 
sition and evaluates the matrix element rea ds 
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